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Abstract. We give the first rigorous (non perturbative) proof of Luttinger liguid behavior in the 
one dimensional Hubbard model, for small repulsive interaction and values of the density differ- 
ent from half filling. The analysis is based on the combination of multiscale analysis with Ward 
identities bases on a hidden and approximate local chiral gauge invariance. No use is done of 
exact solutions or special integrability properties of the Hubbard model, and the results can be in 
fact easily generalized to include non local interactions, magnetic fields or interaction with external 
potentials. 
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1. Introduction 

1.1 Historical remarks. 

The Hubbard model describes electrons in a crystalline lattice, hopping from one site of a lattice 
to another and interacting by a repulsive (Coulomb) force with coupling U > 0. As pointed out by 
Lieb [L] , the Hubbard model in the theory of interacting electrons has the same role of the Ising 
model in the problem of spin-spin correlations, that is it is the simplest model displaying many 
real word features: it is however much more difficult to analyze. It is believed that the Hubbard 
model gives a correct description of the properties of many metals due to the interactions between 
conduction electrons: for instance the phenomenon of Mott transition, the anomalous properties of 
high T c superconductors or the singular properties of quantum wires. However the mathematical 
complexity of the computations is such that this belief is still far from be substantiated by solid 
arguments. While our understanding of the Hubbard model in higher dimensions is really poor, 
the situation is of course better in d = 1; the interest in such a case is not purely academic as 
in this case the model is believed to furnish an accurate description of real systems like quantum 
wires. 

In d = 1 the Hubbard model can be solved exactly by an extension of the Bethe ansatz, as it was 
shown by Lieb and Wu [LW] in the sixties; it was found that in the half filled band case a gap 
opens in the spectrum, so that the system is an insulating, while in the other cases there is no gap 
and the system is a metal. Subsequently many thcrmodinamical quantities were obtained, see for 
instance [O] , [T] ; for a review of exact results see [L] . 

The exact solution is however essentially of no utility for computing the correlations, which are 
the quantities more directly related to physical observables; even if one has the full form of the 
wave functions (what is actually not the case, as the Bethe ansatz gives them as the solutions of 
complicate integral equations), computing the correlations from them is essentially impossible. In 
particular an important question which cannot be answered by the exact solution is if the Hubbard 
model is a Fermi liguid or a Luttinger liquid. The notion of Luttinger liquid was introduced by 
Haldane [H] in the early eighties. While a Fermi liquid is an interacting fermionic system whose low 
energy behavior is close to the one of the free Fermi gas, a Luttinger liquid behaves as the Luttinger 
model; a model describing spinless fermions in the continuum with linear dispersion relation and 
short-range (non local) interaction. The linear dispersion relation has the effect that infinitely many 
unphysical fermions must be introduced to fill the "Dirac sea" of states with negative energy. This 



makes the model a bit unrealistic and of no direct applicability to solid state physics but ,on the 
other hand, the choice of a linear dispersion relation has the effect that, contrary to all other 
models of interacting fermions, the Luttingcr model correlations can be explicitly computed, see 
[ML] . The popularity of the Luttingcr liquid notion increased greatly after the Anderson proposal 
[A] that the high-T c superconductors are, in their normal phase, Luttinger liquids; this proposal 
was based on the conjecture that the Hubbard model in one or two dimensions have a somewhat 
similar behavior, and in particular that they both show Luttinger liquid behavior at least for some 
range of the parameters. Up to now there is no agreement even at an heuristic level on theoretical 
evidence of Luttinger liquid behavior in the d = 2 Hubbard model. On the other hand the Anderson 
proposal stimulated a number of studied about the Luttinger liquid behavior in d = 1 , as a natural 
prerequisite to understand the same question in d = 2. 

Numerical simulations of the correlation functions gave evidence [OS] that the d — 1 Hubbard 
model is indeed a Luttinger liquid; subsequent analytic (but heuristic) results by [PS] and [FK] 
confirmed this result for large U, finding also that the correlations verify an important property, 
the spin-charge separation. For small U the evidence for Luttinger liquid behavior is based on the 
two following facts: 

1) the d = 1 Hubbard model should be equivalent, as far as low energy property are considered, 
to a generalization of the Luttinger model to spinning fermions, the so called g-ology model. 

2) Contrary to the Luttinger model, even the g-ology model is not solvable. However Solyom [S] 
by Renormalization Group (RG) analysis truncated at two loops showed that the g — ology model 
scales iterating the RG toward the Mattis model, a model which is indeed exactly solvable and it 
has Luttinger liquid behavior. 

Given the above two facts, the formulas for the correlations of the d = 1 Hubbard model are 
usually approximated with the formulas for the correlations of the Mattis model, see for instance 
[S]; this is however quite unsatisfactory for a number of reasons: 

1) this approximation means that the effects of the lattice and the corresponding Umklapp scat- 
tering terms in the Hubbard model are completely neglected. There are however strong indications 
that this approximation gives completely wrong predictions at least for properties like the thermal 
or electric conductivity [RA]. 

2) The conclusion of [S] that the Hubbard model scales iterating the RG toward the Mattis model 
is based on a number of peculiar cancellations in the perturbative expansion, checked up to two 
loops. Of course, without an argument stating that the cancellations are present at any order, 
this conclusion is only perturbative and not very solid; if at higher orders the cancellations were 
not present, the effective coupling constants could increase without limit making the analysis 
meaningless. 

As a conclusion, the enormous number of results on the d = 1 Hubbard model can be roughly 
divided in two main classes. A first class are the exact results, based on the Bethe ansatz approach. 
They are (almost) rigorous but they give essentially no informations on the behavior of the corre- 
lations. Moreover they are not very robust, as they rely on delicate integrability properties of the 
Hubbard model and cannot be used to face with apparently harmless modifications of the Hubbard 
model (for instance considering nonlocal but short ranged interactions). The second class of results 
is obtained by a combination of techniques (numerical simulations, bosonization, Renormalization 
Group) and indeed informations on the correlations are found, but the results are not rigorous. 

1.2 The Luttinger liquid construction 

In a series of paper [BM1], [BM2], [BM3], [BM4] a general proof of Luttinger liquid behavior for 
spinless interacting fermions (without any use of exact solutions) has been completed. The con- 
clusion is that interacting spinless fermions are generically Luttinger liquids (independently from 
the dispersion relation, the presence of a lattice, the sign of the interaction etc). A perturbation 
theory for the correlation functions is constructed, and such expansion is not in the strength of 
the interaction but in terms of a set of parameters called running coupling constants, describing 
the effective interaction at a certain momentum scale; the expansion is proved to be convergent 
(and analytic) if the running coupling constants are small, sec [BM1], as a consequence of suitable 
determinant bounds for the fermionic truncated expectations. On the other hand the property that 
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the running coupling constants remain in the convergence radius of the expansion is not trivial at 
all and is due to remarkable cancellations at any order in the expansion. More exactly, the running 
coupling constants verify a set of recursive equations, whose l.h.s. is called Beta function, and their 
boundedness is consequence of dramatic cancellations happening at any order in the Beta function. 
In order to prove such cancellations one decomposes the Beta function in the sum of two terms; 
one, called dominant part, which is common to all spinless d = 1 Fermi system, and the second 
part which depends on the specific model one is considering and which gives a bounded flow once 
one has proved that the dominant part is asymptotically vanishing. The problem is then reduced 
to the vanishing of the dominant part of the Beta function, which coincides with the complete Beta 
function of a model, called reference model, describing interacting spinless fermions with an ultra- 
violet and infrared cutoff. The proof of vanishing of the Beta function is the reduced to the proof 
of suitable (highly non trivial) Ward Identities between the correlation functions of the reference 
model [BM3] , for any value of the infrared cutoff. The problem of implementing Ward identities 
in a model with cutoffs (and in a Rcnormalization Group scheme) is a well known problem in 
Quantum field theory or condensed matter physics. In [BM2], [BM4] a solution for this problem 
was given by finding suitable Correction Identities relating the corrections to the Ward Identities 
due to cutoffs to the correlations themselves. By combining the Ward and the Correction Identities 
the vanishing of the reference model Beta function is proved and the rigorous construction of the 
correlation functions for spinless Luttingcr liquid is then completed. 

Aim of this paper is the extension of the Luttinger liquid construction to the d = 1 Hubbard 
model; such extension is not straightforward as the Luttingcr liquid is not the generic state of 
spinning interacting fermions. The conditions of repulsive interaction U > and not half filling 
must be imposed; technically this is reflected from the fact that the expansions we find cannot 
be analytic in a circle around U = 0, as it would be in the spinless case, and U must be chosen 
smaller and smaller as we are closer and closer to half-filling. We will define an expansion for the 
correlations in terms of running coupling constants, but the presence of the spin increases greatly 
their number; crucial symmetry considerations (based on the SU(2) spin invariance of the Hubbard 
model) and geometrical constraints reduce the number of the effective interactions (which is only 
one in the spinless case) to three (in the non half filled band case) or four (in the half filled band 
case). Again the running coupling constants verify a recursive relation, whose r.h.s. is called Beta 
function, and the expansion is meaningful only if the running coupling constant are small at any 
momentum scale. One can decompose the Beta function in a dominant part and a rest; it turns 
out however that the dominant part is not vanishing. Calling the three (in the not half filled case) 
effective interactions g\, gi, 54, it turns out that, truncating the beta function at the second order, 
gi tends to vanish (if U > 0) while <72,</4 remains close to their initial value. In order to prove 
that such a result is valid non perturbatively, that is including all orders contributions, one has 
to prove a property which we will call partial vanishing of the Beta function. Such property is 
derived by a suitable reference model, which verifies formally (if cutoffs are neglected) proper gauge 
symmetries. Quite surprisingly, the cancellations on Beta function of the Hubbard model, which 
verifies an SU(2) spin symmetry, will be obtained by a reference model not SU(2) invariant. We 
derive suitable Ward and correction identities for the reference model, and by them the partial 
vanishing of the Hubbard Beta function is proved. Hence the running coupling constants are small 
if Uis small enough and the rigorous construction of the correlation functions for Hubbard model 
is completed. The analysis can be easily generalized to include the presence of a magnetic field or 
non local interactions. 

1.3 The Hubbard model 
The Hubbard model Hamiltonian is given by 

h = -t (°i a+i,<x + a i+i,a a x,a) + u Yl a t.,+ a x,+ai,-a~,- - it Y a t,* a ~x,o (1-1) 

xeA&=± xeA xeA&=± 

where A is an interval of L points on the one dimensional lattice of step 1, which will be chosen 
equal to (—[L/2], \(L — l)]/2) and a~ a is a set of fcrmionic creation or annihilation operators with 
spin a = ± satisfying periodic boundary conditions; t — 1/2 is the hopping parameter, U > the 
coupling and /i the chemical potential. The Hamiltonian verifyes an SU(2) spin symmetry. 
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A generalization of the above hamiltonian including the effect of a short-ranged (instead of a 
nearest-neighbor) interaction, and the presence of a magnetic field, is the following 

H = ~\ J2( a t° a x+h<J + a t+l.a a x.a) +UJ2 V ( X - V) a i., + a x., + a t^ a y,- ~ V^2 a tv a x,a ( L2 ) 

x,<j x,y x,a 

+h(a+ +a - + - a+_a-_) 

When the interaction is given by Uv(x — y) = U5 x-y + V6 x +i iV the model is known as the U — V 
model. 

We consider the operators = e Hx °a^ a e~ Hx °, x = (x,Xo) and — (3/2 < x n < 0/2 for some 
(3 > 0; on x , which we call the time variable, periodic boundary conditions are imposed. 

Many physical properties of the fermionic system at inverse temperature /? can be expressed 
in terms of the Schwinger functions, that is the truncated expectations in the Grand Canonical 
Ensemble of the time order product of the field at different space-time points. If 

Tre~ l3H X 

< X > L , = Tre _ fjH (1.3) 

the Schwinger functions are defined as, if e = ± 

Sx,/3(xi,£i,0-i;...;xi,£i,<7i) =< a^^.-.a^^ > Ly0 (1.4) 

We will denote by 5(xi, ...) the lirn^^oo of (1.4). An important role is played by the two point 
Schwinger function 

<Sx,/3(x, +, <r; y, -, a) = S LJj (x, y) (1.5) 

Denoting by Sx,/3(fc, #0) the Fourier transform of iSx, )( g(x) with respect to the x variable, = 
SL,p{k, 0~) is the occupation number, the average number of particles with momentum k. Another 
important physical quantity is the density-density correlation function 

^x,y =< PxPy > - < Px >< Py > (1-6) 

1.4 The non interacting U = case. 

The two point Schwinger function in the non interacting case is given by, in the limit L,f3^>oo 

/co pit -ik(x-y) 

dk ° dk -ikn + U-cosk (L?) 
-oo J -n '"O T /I COb fi. 

where k = (k, ko). It is easy to check that one can write, if u. = cosp° F and vo = sinp° F 

p iup° F (x-y) 

S (x,y)=V -— +.9(x,y) (1.8) 

^— ' v x + iwx 

UJ — ± 

with |<?(x, y)| < 1+ | x fy|i+a , a positive constant; that is, the two point Schwinger function decays 
as 0(|x — y| _1 ) oscillating with period Important physical properties are: 

1) the occupation number is given by — x(|k| < p F ), that is is discontinuous. 

2) The bidimensional Fourier transform of density correlation has singularities at (±2p F ,0) and 
(0, 0); in (±2p F , 0) it has a logarithmic singularity while in (0, 0) the Fourier transform is bounded. 

3) The one dimensional Fourier transform at xo = of the density correlation is continuous, while 
its first derivative in k has a first order discontinuity in k = 0, ±2p F . 

1.5 Main result 
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Our result can be informally stated in the following way 



In the not half filled band case and in the weak coupling regime, the (repulsive) Hubbard model 
(1.1) is a Luttinger liquid. 

A more formal statement is the following theorem. 

1.6 Theorem 1. Consider the hamiltonian (1.1) with — 1 < n < 1 and fi ^ (not filled or half 
filled band case); there exists an e > such that, for < U < e 

a) the two point Schwinger function (1.5) is given by, in the limit L,(3 — > oo 

S ^ y) = E - ( — — t I,, + *(x, y) (i.9) 

v(x - yo) + iw(x - y) |x — 

with 

n = all 2 + U 2 f (U) p F =cos-V + /i(£0 v = v + f 2 (U) (1.10) 
where a>0, \fo(U)\, \fi(U)\, \f 2 (U)\ < CU and 

\d^d2A{*,y)\ < CU- l —— |S(x,y)| < , ° .... (1.11) 

/or suitable positive constants C, 6, if d denotes the discrete derivative. Moreover the occupation 
number nk is continuous atk = ±pp but its first derivative diverges atk = ±pp as \k— (±pp)\~ 1+ ' n . 

b) The density- density correlation function (1.6) can be written as 

ft x = cos(2p F x)Q a (x.) + ft b (x) + ft c (x) , (1.12) 

with ^ 4 (x) 7 i — a,b,c, continuous bounded functions, which are infinitely times differ entiable as 
functions of xq, if i = a,b. Moreover 

n-(x)= 1 + ^ l(x) 



2n 2 [x 2 + (vx o n+m ' 
Qb (x) _ 1 ( x 2 -(x/v ) 2 } 

|4(x)| < ^ I^WI < 1 + |x C y|2+e , (1-14) 
for some constant C, where rji = —bll + U fi{U) with b > and \f&{U)\ < CU. Finally for a = a,b 

|^^n»-(x)| < |x _ y| ? + „ 0+ni (i.i5) 

c) Let O(k), k = (k,k ) G [— 7r,7r] x R , ffte Fourier transform o/0(x). For any /ixed k u>zi/i 
k 7^ (0,0), (±2pi?,0), fi(k) is uniformly bounded; moreover 

|A(k)| < c 2 [l + t/log Ikr 1 ] (1.16) 

near k = (0, 0), and, at k = (±2j>f, 0), diverges as 

|k-(±2 PF ,0)| 2 "7lml (1-17) 
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Let G(x) = Q(x,0) and G{k) its Fourier transform. Then dkG(k) has a first order discontinuity 
at k = 0, with a jump equal to 1 + 0(U), and, at k = ±2pp, diverges as \k — (±2pp)\ 2ni . 

Remarks 

a) A naive estimate of e in the above Theorem is e = 0{\(i\ a ) for some constant a, for /j, close to 
0; that is U must be taken smaller and smaller as we are closer and closer to the half filled band 
case. 

b) A first effect of the interaction is that the Fermi momentum pp is modified by the interaction 
by 0(U) terms. 

c) More dramatic is the effect of the interaction on the long distance asymptotic behavior of the 
physical observables; it turns out that the two point Schwinger function decays faster in presence 
of the interaction, while the correlation function decays slower. The large distance decay is power 
law with anomalous critical indexes depending non trivially by the coupling U. 

d) As a consequence the occupation number rife, which in the non interacting case have a discon- 
tinuity at k = ±Pf, has no discontinuity in presence of the interaction; this proves that the d = 1 
Hubbard model is a Luttinger liquid in the sense of [H] . The lack of discontinuity in the occupation 
number can be physically interpreted saying that fermionic quasiparticles are not present. 

The interaction changes the log-singularity at k — (±2pp,0) of the bidimensional Fourier trans- 
form of the density correlation in a power law singularity, with a nonuniversal critical index 0(U). 
This enhancement of the singularity is considered a signal of the tendency of the system to develop 
density wave excitations with period tt/pf, generically incommensurate with the lattice. On the 
other hand the singularity in k = (0, 0) is much weaker, that is at most logarithmic. 

In the same way, the interaction leaves invariant the singularity of the first derivative of the 
one dimensional Fourier transform of the correlations in k = (a first order discontinuity) while 
the singularity in k = ±2p F is changed by the interaction from a discontinuity to a power law 
singularity. 

e) The two points Schwinger function and the density correlation can be written as sum of two 
terms; one which is very similar to corresponding quantities in the Luttinger model, and in which 
the dependence from pp is quite simple (they can be written as oscillating terms times terms 
which are free of oscillations, in the sense that each derivative increases the decay by a unit, see 
(1.11), (1.15)) and another (non Luttinger like) in which the dependence on pp and the lattice steps 
is very complicate; this last term decays faster than the Luttinger like terms but the derivatives 
do not increase the decay for the presence of oscillating terms. The non Luttinger like terms have 
Fourier transform which is bounded; however sufficiently high derivatives of the Fourier transform 
can be singular for values different respect to k = 0,±pf,±2pf (such singularities were indeed 
observed in numerical simulations, see [OS]). 

d) Our results provide a proof of Luttinger liquid behavior, but they are still not enough ac- 
curate to prove an important property called spin-charge separation, which is believed true for 
the Hubbard model; namely that the asymptotic behavior of the two point Schwinger function is 
(x + iv c x)~ 2~ Vc (xo + iv s x)~^~ T,s , with v c — v s = 0{U) and i] c ,r/ s = 0(U); (1.9), (1.11) is compat- 
ible with such behavior but is not enough accurate to prove it. Another property which could be 
probably proved by an extension of our techniques is the Borel summability of our critical indexes 
as a function of U. 

e) Finally we could consider a short range instead of local potential, that is (1.2) with h = 0. In 
such a case the condition U > is replaced by the condition Uv{2pp) + F{U) > 0, where F(U) 
is a suitable 0(U 2 ) function. Note that the linear term is vanishing for sufficiently long range 
interactions such that v{2p F ) = 0. 

1.7 The Hubbard model in a magnetic field 

Let us consider the Hamiltonian (1.2) with h ^ 0; the presence of a magnetic field destroys the spin 
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rotation invariance. Moreover it turns out that that one can consider also attractive interactions, 
if the interaction is smaller than the magnetic field. Calling S ai L,p{x, y) =< Vx.tr^y.a >l,0 we 
prove the following result. 

1.8 Theorem 2. Consider the hamiltonian (1.2) with —1 < ji < 1 and < h < ho for a suitable 
constant ho ; assume also that cos (/i + h) + cos (/i — h) ^ n. There exists positive constants 
£\,S2 ( depending on p and h,and £2 vanishing as h — > 0) such that, if —£2 <U<£\ the two point 
Schwinger function is given by, in the limit L, (3 — > 00 

S CT (x,y) = ]T- . r — |- 1 ' + ^(x, y) (1.18) 

^w(a;o-2/o)+«w(a;-y) |x - y^ 



r? = a£/ 2 +0([/ 3 ) p£ = cos-V+sign(cr)/i)+0(LT) i> = sin(cos _1 (At+sign (a)h))+0{U) (1.19) 
where a > and 

|cWM(x,y)| <C£7|x-y|- n "-^ |S(x,y)| < 1 + |x C y|1+9 (1.20) 

/or suitable positive constants C,9. 

The other statements in the previous Theorem can be repeated with some obvious modifications. 
The above result says that the Hubbard model is still a Luttinger liquid even in presence of 
a magnetic field; this happens even in the attractive case, if the interaction is smaller than the 
magnetic field. 

1.9 Contents 

In §2 and §3 we write the Hubbard model (1.1) partition function as a Grassmann integral, and 
we define a multiscale integration procedure; we get an expansion in terms of running coupling 
constants, whose regularity properties are stated in Theorem 3. In §4 we study the flow of the 
running coupling constants and in §5 and §6 we derive the cancellations of the Hubbard model 
Beta function by Ward identities and Correction identities of a suitable reference model. Finally 
such results are applied in §7 to the computation of the Schwinger functions and the correlations 
and in §8 the presence of the magnetic field is included. We rely on many technical results already 
obtained in [BM1-4] (the presence of spin has a small effect on the proof of convergence, for 
instance) and we focus mainly on the difference with respect to the spinless case. 



2. The ultraviolet integration 



We assume \x € SI n (—1, 1), where VL C — {pL : \fi ± 1| < fi, where /x > is a fixed constant. 
We call V LJ3 = V L x V p , with V L = {k = 2?rn/i,n e Z,-[L/2] < n < [(L — l)/2]} and 
V f3 = {ko = 2(n + l/2)7r//3, ne Z, -M < n < M - 1}; moreover we define 

1=1 — 5 tcospF = [i — v (2.1) 



with 5, v suitable counterterms to be fixed properly in the following. We introduce a finite set 
of Grassmanian variables {-0^ }, one for each k € T^L,p-, and a linear functional P{dijj) on the 
generated Grassmannian algebra, such that 



/ 



Wfe^k^ = £/3<5k llk2 <W 2 5(ki) , m = - —~ 1 ~ ■ (2.2) 

1,1 -ik + tcospF — t COSK 



We will call g(k) the propagator of the field. 
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We define also Grassmanian field ip^ i s defined by 

<° = h S ( 2 - 3 ) 

such that 

1 £ e -^-y)j( k ) = / F ( W + E ^( x;y)i (2.4) 

It is well known that the partition function Z =< e~@ H > L,p can be rewritten as the limit M — > oo 
of the Grassmann integral 

y P(dtP)e~ v (2.5) 
where P(dtp) is the Grassmann integration with propagator (2.4) and 



dx Q V'i+C+V'x.-C- + " / ^oV V4.<xC<x + 5 / dx J! k.B^U^x.a 
■PI* x J -PI2 x . <7 J -P/2 x.y,* 

(2.6) 

where t x>y = \5 y ^ x+ \ + j5 x , y +i and P(dip) has propagator given by <?(x, y) (2.4). Let T 1 be the 
one dimensional torus, ||fc — fc'|| T i the usual distance between k and k' in T 1 and ||fc|| = ||fc — 0||. 
We introduce a scaling parameter 7 > 1 and a positive function x(k') <E C 00 ^ 1 x R), k' = (fc', fco), 
such that 

X (k')=x(-k') = (i |f|k;i<*o^aoV7, (27) 
AV ; AV ; [0 if |k'| > a v Q , v ; 

where 

\V\ = ^/k 2 + (v \\k'\\ Tl ) 2 , (2.8) 
a =min{p F /2,(7r-p F )/2} , (2.9) 

The definition (2.7) is such that the supports of x(fc — pf, ko) and \(k + Pf, ko) are disjoint and 
the C°° function on T 1 x R 

U.v. (k) = 1 - X (k - PF, k ) - x(k + Pf, k ) (2.10) 

is equal to 0, if [i>o||(|fc| - Pf)\\tA 2 + H < t%. We define 

ff L ^(x;y)= 3 "-(x,y)+ 3 "-(x,y) (2.11) 

with 

9 V 



ip, ' — ik — t cos k -\- 1 cos pf 



^. (X;y) 1 E e -, k(x - y) n r±1 x(fc-^,fa) 

Lp i —f —ik — tcosk + tcospF 



(2.12) 



From the integration over ^("-"-J we get 

e -^,, = g-^^+to y p( # (i.r.)) e -v'°>Wf'-'-») 7 V (0) (0) = ^ (2-13) 



e 



P(#("-"-))e- v ^ u r +^ ; ) . (2.14) 
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We will call = V> ( - 0) ; V( )(^ ( - 0) ) can be written in the form 



In 



v (0 W* 0) ) = E(zkE E IlCf*- 

n=l y t ' a ki,...,k 2 „ i=l 
2n 

• ^i^Jk^.^k^O^E^k,), 

i=l 



(2.15) 



where a = (a\, . . . , <T2„), w = (u>i, . . . , W2n) and we used the notation 

5Ql) = 6(k)6(ko) , 6(k) = Lj2h,2nn, S{ko) = p6 kofi • (2.16) 

We prove in the Appendix that 

IW^^ki, ...,k 2 „_ 1 )| < L/3Cmax([/, \ V \)™*V.»M (2.17) 
The SU (2) spin invariance implies that the quartic terms have the following form 

77^4 E W 4 ( °i(ki,---,k4)E 5 (E £ * ki )[^ 

^ ' ki...,k 4 (T i 

(2.18) 

where is spin independent. 

3. The infrared integration 



3.1 Quasiparticles. 
We define also, for any integer h < 0, 

A(k') = x(7^k') - x (7- fc+1 k') ; (3.1) 

we have, for any ft, < 0, 

X (k')= E hW+xh-^'). (3.2) 

h=h+l 

Note that, if h < 0, / h (k') = for |k'| < toT^ 1 or |k'| > t 0l h+1 , and / h (k') = 1, if |k'| = t j h . 
We finally define, for any h < 0: 

fh(k) = fh(k-pF,k ) + f h (k+p F ,k ) ; (3.3) 

This definition implies that, if h < 0, the support of /h(k) is the union of two disjoint sets, and 
A^. In A^, fc is strictly positive and ||fc — pfWt 1 < a o7 h < a o, while, in , fc is strictly negative 
and ||fc + pfIIt 1 < &07 ■ The label /i is called the scafe or frequency label. Note that, if k € 2?l,/3, 
then |k± {p F ,0)\ > + (vqttL- 1 ) 2 , by the definition of V L . fi . Therefore 



A(k) - Vh < h LtP = mm{h : t Ql h+1 > Vt^ 1 ) 2 + {v^L^f} , (3.4) 
and, if k £ T^l,0, the definitions (2.10) and (3.3), together with the identity (3.2), imply that 

o 

1= E A (k) + /;.,. (k) . (3.5) 

K—Klq 



2/ febbraio/ 2008; 13:13 



9 



We now introduce, for any h < 0, a set of Grassmann variables ip^, u such that 

J P(d^)^^+W ^ = LflS^S^d^g^Hk',) . (3.6) 

where 

g W(k / )= - Ih^lM (3.7) 

' w — iko — tcos(fc' + wpi?) + t cospi? 

We introduce also the Grassmann variables 

tf$& = T5 £ ^t e±lk '' x ( 3 - 8 ) 

so that 

' = <W Ws^fry) • (3.9) 



where 

ff W(x,y) = _L ^ e - k '( x -y) ff W(k') • (3.10) 

It holds that 

fp(#("'»)cf-)^ ( ;;^v E E^^'^W)- ( 3 - u ) 

h=h L ,f) w=± 

The above identity implies that, if F(tp^- r ^) is any analytic function of the variables ^ l - r ) 

/ P(#("V» (ir) ) = / f[ / PW W )F( E E e "" PFX ^) (3.12) 
We define also 

ft 

C,- 1 (k)= E AW (3-13) 

fc= — oo 

3.2 Multiscale integration. The integration of the infrared part is done in an iterative way. 
Assume that we have integrated the scales 0,— 1, ft + 1 and that we have found 

J Pz h ,c h (d^- h) )e- V<hH ^ ( - h)) - LflEh , V< h >(0) - , (3.14) 

where 

p Zh ,cM^ h) )=K- 1 n n <fw k -g h) 

k'enw=±i 

exp-^ E ^WSt-^ + ^^^ + ^Mcos^-l))^ 1 (3-15) 

C- 1 (k)>0 

and 

In 



^(^^-Etz^e e ncir 

• MnljK, -Xn-l) ^E ^ + E e * W ***) 



^ (L/3) 2 " ^ ^ 11 

- k i-- k »»* (3.16) 

2n 2n 



»=1 i=l 
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and in particular the quartic terms have the following form 

E E<2( k i'-' k 3)E^(E £ ^+E £ ™) 

^ ' ki,..,k^ « <r i i 

[V'k 1 , Wl ,crV'k2,a)2,o-^'k3,W3,o-V'k4,W4,cr + ^ki ,wi ,cr V'ka ,w 2 ,<r ^k 3 ,w 3 V'ki ,W4 (3-17) 

Note that there exists a scale h such that, for h < h are present in (3.17) only the monomials 
verifying 

4 

HE e < w <PHIn =0. (3.18) 

i=l 

In fact by the compact support properties of the propagators ll^^i^illT 1 < 4ao«o7 /l+1 an d if 
(3.18) is not satisfied 1 1 X)i=i ^^P-H It 1 > C|pf — f |; hence h — 0(log\pp — § |). The condition 
> surely implies that |pf — f | > (hence that h is finite if \i ^ 0), for {/ small enough (than 
O(p)), as = cos -1 fi + 0(U). 

3.3 TTie localization operator. 

We split the effective potential as jCV^ + TZV^ h \ where 1Z = 1 — C and C, the localization 
operator, is a linear operator on functions of the form (3.16), defined in the following way by its 
action on the kernels W^nu- 

1) If 2n = 4 we define 

£<L( k 'i' k 2> k 3) = L- 1 8(^e I Lo lPF )wl>ljk ++ ,k ++ ,k ++ ) , (3.19) 

»=i 

where k,,,,/ = (tj-kL^ 1 , r/'irfl^ 1 ). 

2) If 2n = 2 (in this case there is a non zero contribution only if u)\ — L02) 

CW£IM = \ E ^2°L(kw){l + ^+V^o| , (3.20) 



3) In all the other cases 



^l,o,(ki,.-.,k / 2 „_ 1 )=0. (3.21) 



In the not half filled band case pp 7^ § the condition 5(^ i=1 e^ipp) ^ is equivalent to the 

condition X^i=i eiCLli 7^ 0- Then the action of C if n = 2 is non trivial only if Ylt=i £ i LjJ i = an d 
there arc only the following possibilities for u>\, u>2, W3, W4: 

— w, — w); (w, — w, — u),u); (ui, ui, w, u) (3.22) 

In the half filled band case pp = j the action of £ is non trivial also if loi = U3 = —lo 2 = —0J4. 
We get 

rV< h >(V0=7 fc nfc^ fc ^)+*fc^^ (3-23) 
where 



= W E E ^k',os^k',osa 



k' OJ.CT 
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Fa = 4lY^^ SinPFS[nk/ + COSp F (cOsfc' - l)]^^k',a,,^k',a;,a 
k' u,a 

1 ' k'j.k^k^k^ w ct.ct' » 

^ ' k'j.k^k^k^ u,o,cr' i 

= - — — — 7 V V(5(V £j k')?A+ Vr7' >Vv - 

Note that 

72, h = {u>PF, up F , -U>p F , -UPf) 7l,h = (<^PF, -WPf, -WPf, ^Pf) 

74, h = W w (ljpf^Pf,^Pf,^Pf) (3.24) 

and in particular 

74 :0 - Uv(0) + 0(U 2 ) 72,0 - Uv(0) + 0(U 2 ) 71,0 - Uv(2 PF ) + 0(U 2 ) 

In the case of local interactions v(p) = 1. Note also that the spin symmetric part of g 4 is vanishing 
by Pauli principle. 



3.4 Renormalization. 
We write (3.14) as 

J P z ^ C(i (^(<^)) e -^< h »(^<^»)-KV^(^ ( ^»)^L/3 Sh f (3 25) 

and we include the quadratic part of LV^ given by J dk' ^ ^ u a (—iko+uj sin fc'+cospir(cos fe'- 
i))^k' w CT in the free integration; we call 

CV h = CV {h) - z h j dk'J2^^A- ik o + ^ s ^ k ' + cos PF(cosk'-l))^ loja (3.26) 

J UJ,<7 

so that we obtain 

/ P^_ i Ch (#(^ fe ))e-^' l (^ ( -' ,) )- w( '' > (^ ( ^))-i/3i5 ) , 7 (3.27) 

where 

Zfc-i(k) ='^(1 + ZhCft *(]£)) (3.28) 
It is convenient to rescale the fields: 

+^F a ( v /^I 1 > ( - /0 +7V^(V^T* ( - h) ) + ftV< h >(V^hV' ( -' ) , (3.29) 



where 

f/i = TT~~ n h S h = — [an - z h ] g i>h = [— -} 2 -fi,h (3.30) 

Ah-l ^h-l ^h - 1 
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Finally the r.h.s. of (3.27) can be rewritten as 

e"^ / Pz^c^W^*) J F.^j-.^le 41 "^^ 1 ) , (3.31) 



where 



and 



Z h -i = Z h {\ + z h ) /,(k') = A(k')[l + ^^y] (3.32) 



f P - i(dibWU~W ib +{h) -5 5 ^ l \^,y) 



1 1 _ — . e -»k'(x-y) 

/ / (^) 33) 



Zu_\L8 ^ ' — ik + u)v sinfc' + cosp^fcos k' — 1) 
We then integrate tp^ 

J = e - v<h " 1) (V^^ < - h) ) (3.34) 

and the procedure can be iterated. 

Note that the quartic terms in CV h can be written in coordinate representation in the following 

way 

J dx 9%[i>tu ) ,*i J x,u>,ai J t-^.,^x,-^,v + ^ 



LU,(T 



+9l K^.^x-^^x-^^x.b,,, + V^,^x ,- W ,a^x -a^x.w ( 3 - 35 ) 
+94[^x", w ,<t^x 2 , w ,,t^x, w -a^x, w - J 

where J g?x = J dxo J2 X ■ Finally note that the propagator is written as 

ff£(x - y) = fli fc ,i(x - y) + rf, L (x - y) (3.36) 

where 



1 p -ik(x-y) 

^ i( x-y) = — ^/,(k)- 



/3L ^ — ik + Luk 



and for any positive integer N 



\r h {x-y)\<C NZ 



l + ^lx-yl)^ 



It is easy to verify that gj/ verifies the same bound of r/j with a "f h less. We call v k — (ffc, 5k, gi,k, 92, k, 54. k), 
k < and vi = (v,S,U); moreover we call g k = (gi,k, g2,k94,h) and fj, k = (52,fc,ff4,fe), fc < 0. The 
above integration procedure generates a power series expansion for w( x i' ••' x 2n) in (3.16) in 

terms of the running coupling constants Vk, k = 1, 0, — 1, — 2, h, which is indeed convergent if 
they are small enough. More exactly it holds the following result. 



3.5 It holds the following crucial result. 

Theorem 3. Assume that fj, ^ 0, ±1 and sup fc>?l \vk\ < £h! assume also that, for some constant 
c, sup fc>/l — e<:e S' then there exists e such that, for Sh < £ the functions W^^^xi, ...x„) are 

analytic in the running coupling constants (v k )k>h and, for a suitable constant C,a 

J d Xl ..d Xn \W^,J^u ...x„)| < (Ce h h a r^ n ^Lp^ 2 -^ h (3.37) 
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Sketch of the proof. The proof is essentially identical to the one of Theorem (3.12) of [BM1] about 
the spinless case. The only important difference is that there exists a finite scale h = 0(log \pf~ § |) 
such that for h < h there are no contributions to the effective potential V h (3.16) with n = 2 and a 
choice of uj,e such that (3.18) is not verified. We can write the effective potential V^ h \\ZZh^- h ^), 
for h < 0, in terms of a tree expansion, similar to that described in [BM1]. 




Fig. 1 



3 j + 1 K -1 +1 

Fig. 1: A tree r and its labels. 



We need some definitions and notations. 

1) Let us consider the family of all trees which can be constructed by joining a point r, the root, 
with an ordered set of n > 1 points, the endpoints of the unlabeled tree (see Fig. 1), so that r is 
not a branching point, n will be called the order of the unlabeled tree and the branching points 
will be called the non trivial vertices. The unlabeled trees are partially ordered from the root to 
the endpoints in the natural way; we shall use the symbol < to denote the partial order. 

Two unlabeled trees are identified if they can be superposed by a suitable continuous deformation, 
so that the endpoints with the same index coincide. It is then easy to see that the number of 
unlabeled trees with n end-points is bounded by 4™. 

We shall consider also the labeled trees (to be called simply trees in the following) ; they are defined 
by associating some labels with the unlabeled trees, as explained in the following items. 

2) We associate a label h < with the root and we denote Th,n the corresponding set of labeled 
trees with n endpoints. Moreover, we introduce a family of vertical lines, labeled by an an integer 
taking values in [h, 2], and we represent any tree r € Th,n so that, if v is an endpoint or a non 
trivial vertex, it is contained in a vertical line with index h v > h, to be called the scale of v, while 
the root is on the line with index h. There is the constraint that, if v is an endpoint, h v > h + 1. 

The tree will intersect in general the vertical lines in set of points different from the root, the 
endpoints and the non trivial vertices; these points will be called trivial vertices. The set of the 
vertices of t will be the union of the endpoints, the trivial vertices and the non trivial vertices. 
Note that, if v\ and V2 are two vertices and Vi < v%, then h Vl < h V2 . Given a vertex v, which 
is not an endpoint, x„ will denote the family of all space-time points associated with one of the 
endpoints following v. Moreover, there is only one vertex immediately following the root, which 
will be denoted vq and can not be an endpoint; its scale is h + 1. Finally, if there is only one 
endpoint, its scale must be equal to +2 or h + 2. 

3) With each endpoint v of scale h v — +2 we associate one of the three contributions to V(V^-^), 
written as in (2.6) and a set x„ of space-time points, the corresponding integration variables. With 
each endpoint v of scale h v < 1 we associate one of local terms in LV lyhv ^ 1 ^ (3.29); we will say that 
the endpoint is of type g\, gi and so on depending on the term we associate to it. 

Moreover, we impose the constraint that, if v is an endpoint and x„ is a single space-time point 
(that is the corresponding term is local), h v = h v > + 1, if v' is the non trivial vertex immediately 
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preceding v. 

4) If v is not an cndpoint, the cluster L v with frequency h v is the set of endpoints following the 
vertex v; if v is an endpoint, it is itself a (trivial) cluster. The tree provides an organization of 
endpoints into a hierarchy of clusters. 

5) We introduce a field label / to distinguish the field variables appearing in the terms associated 
with the endpoints as in item 3); the set of field labels associated with the endpoint v will be called 
I v . Analogously if v is not an endpoint, we shall call I v the set of field labels associated with 
the endpoints following the vertex v; x(/), <r(/) and to(f) will denote the space-time point, the a 
index and the u) index, respectively, of the field variable with label /. 

If h < 0, the effective potential can be written in the following waysee [BM1]: 

oo 

V^(^^)+LpE h+1 = J2 E VW(r,y/Z^^ h) ) (3.38), 

n=l reT hj „ 

where, if v is the first vertex of t and ti,..,t s (s = s Vo ) are the subtrees of r with root vq, 
V( h \r, y/Zhipt-' 1 ^) is defined inductively by the relation 

v^(T,^z~ h ^ h y) = 

IZ^^^+I)^ ^(<W-i) ); .,V^(r s , y/Z^W)] , (3 ' 39) 
and F( fc + 1 )(T i ,V^ ( - fc+1) ) 

a) is equal to KV (h+1 *>(n, y/ r Z^ip ( ^ h+1 '>) if the subtree n is not trivial; 

b) if n is trivial and h < — 1, it is equal to one of the terms in LV^ h+1 ^ (3.29)or, if h = 0, to one 
of the terms contributing to V(V'- 1 ) (2.6). 

It is then easy to get, by iteration of the previous procedure, a simple expression for V^ h \r, y/Zhipt-^), 
for any r e %i, n - 

We associate with any vertex v of the tree a subset P v of I v , the external fields of v. These subsets 
must satisfy various constraints. First of all, if v is not an endpoint and v\, . . . , v Sv are the vertices 
immediately following it, then P v C UiP Vi ; if v is an endpoint, P v — I v . We shall denote Q Vi the 
intersection of P v and P Vi ; this definition implies that P v — DiQ Vi . The subsets P Vi \Qvi, whose 
union will be made, by definition, of the internal fields of v, have to be non empty, if s v > 1. 

Given r G Th. n , there are many possible choices of the subsets P Vl v € t, compatible with all the 
constraints; we shall denote V T the family of all these choices and P the elements of V T . Then we 
can write 

V (h \r,^Z~ h ^) - E ^ ( ' l) (^P) 5 (3-40) 

Calling W^p the kernels of V^ h \r, P) (see (3.16)) and repeating the analysis in §3 of [BM1] one 
gets the following bound (analogous to (3.105) of [BM1]) 



dx„\W. < C"LBely- hD '' p -> ■ 

tjnot e.p. v ' 

where z(P„) = 2 if |P„| = 2 and z(P v ) = 1 if \P V \ = 1 and || X)/eP v e (/) w (/)?'i :, l It 1 = 0; moreover 
x(Pv) are defined so that x(Pv) = if \P V \ = 4 , h v <h and 1 1 X)/eP„ £ (/) w (/)Pf| It 1 7^ 0, and 
x(P„) = 1 otherwise. 

We call 2 — — z(P„) the dimension of the vertex w in the tree. If no renormalization is defined 
1Z = 1 then one gets a similar bound with z v (P v ) = 0. Hence if 1Z = 1 the vertices w with \P V \ =4 
have vanishing dimension (marginal terms) while if \P V \ = 2 they have positive dimension (relevant 
terms). The presence of the \~ functions in (3.41) is easily understood by noting that one can insert 
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freely such % functi ons in momentum space, then one passes to coordinate space and make bounds 
using the Gram-Hadamard inequality as in [BM1]. 

For any v such that h v < h it holds — 2 + ^f^- + z(P v ) > 1, that is the dimension is negative, while 
if h v > h it holds -2 + ^4 + z(P v ) > 0. 

We have to perform the sums over r and P. The number of unlabeled trees is < 4"; fixed an 
unlabeled tree, the number of terms in the sum over the various labels of the tree is bounded by 
C™, except the sums over the scale labels and the sets P. 

In order to bound the sums over the scale labels and P we first use the inequality, for a constant 
< c < 1 



n x(Pv)(z hv /z hv _ 1 ) 



\Pv\/2 _ r _ OJ . JSl 



7 



[-2+^+z(P v )\ < 



-not e.p. (3 42) 

<[]l(x(hv<hh-^- h ^+ X (h,>h))}[ II X(l^l>4)7"^], 

v v not e.p. 

where v are the non trivial vertices, and v' is the non trivial vertex immediately preceding v or the 
root. Then it holds that, noting the the number of nontrivial vertices is bounded by n 

J2 ll[(x(hv < hh~ c ^- h ^ + X (hy > h))} < C n \h\ an (3-43) 

{ha} v 

for some numerical constant a. Finally the sum over P can be done as described in [BM1]. ■ 
Remark 

By (3.42) we get also that the bound for a tree r € T h ^ n with at least a vertex at scale k improves 
by a factor <-y e (' l - fe ); this property is called short memory property. 

4. The flow equation 



4.1 Second order analysis. 

By the iterative integration procedure seen in the previous section it follows that the running 
coupling constants verify a recursive relation whose r.h.s. is called Beta function: 



Z h - 



- = 1 + z h {v h ,..,vx) (4.1) 



Z h 

Vh-i=lVh + f3i h) {v h ,.., Vl ) (4.2) 

5 h -i = 5 h + pf\v h ,..,v 1 ) (4.3) 

9i,h-i=9i,h+^ g h h v h,.;Vi) (4.4) 

with i = (1, 2, 3). The above equations are also called flow equations. The functions Zh, Pi h \ 0^ P^}) 

are expressed by the tree expansion seen in §3 (for details, see [BM1]). The contribution to f3^\ 
from the trees with two end-points associated to the quartic running coupling constants is given 
by, if J* dr = J^ 2 /2 dr 

^4 [drgl k \r)gW(-r)g hk g hk =4 f dv £ g^ {v)g^l{-v)g hh g hk (4.5) 

k<h •* J k=h,h+l 
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Using that g\ t h — #1,^+1 = 0(v%) and (2.11), and computing the equation analogue to (4.5) for 
g2,h-i an d 9A,h-\ we get that g^h verify the following equations 

91, h-i = gi,h - agl h + 0{v 2 hl eh ) + 0{vl) 

92, h-i = 92,h ~ \9lh + O(vh 0h ) + 0(4) (4.6) 

94,h-i =g 4 , h + 0(v 2 hl eh ) + 0(vl) 
with a a positive constant, given by 

a = a 1 + a 2 = 4j dr[g^(r)g^Jr) + ^(r)^(r)] (4.7) 

If we neglect the cubic contributions O(v^) it is easy to see that the flow is bounded (in sense that 
the quartic running coupling remain smaller than 0(U) for any h) if U > 0; in the general case 
in which the interaction is non local the conditions is 51,0 = Uv(2pp) + 0(U 2 ) > 0. Of taking 
into account all higher order terms could destroy such behavior; aim of the following sections is 
to prove that also taking into account the full Beta function the quartic running coupling remain 
smaller than 0(U). 

4.2 Beta function decomposition. 

We have two free parameters pot our disposal, v and <5; we will show that we can fix them so 
that v h = 0(U 2 -f Th ) and S h = 0(U 2 -f Th ). We fix then our attention on the flow equation for 

9l,h> 92,h,94,h- 

More explicitly (4.4) can be written as 

9im-i = 9i,h - gi,h[ a i9i,h + a 2 gi,h+i] + G\(g h , .., g ) + ^ g 1 , k gi,k' H h,k,k' ( v h, ■■, v ) + Rl(vh, -~vi) 

k,k' 

92M-1 = g2j l -^gi,h[aigi,h+a2gi,h+i]+(3h(t J -h, ■■^o)+G 2 h {g h , .., go)+^gi, k gi, k > Hi kk ,(v h , ..,v )+R 2 h (v h , ...v^ 

k,k' 

(4-8) 

9i.h-\ = 92,h + Phinh, -,Vo) + Gl(g h , ..,g ) + ^2gi,kgi,k> H h,k,k>( v h, -,v ) + Rt(v h , ...vx) 

k,k' 

where the following definitions are used: 

1) We write in (4.1) z k — z\ + z 2 , where z\ is defined iteratively as the sum of all trees with 
only end-points at scale < and with propagators g^\ see (2.11), and in which , k' > k is 
replaced by 1 + z\, . 

2) The functions f3^, (3^, G^, G\, G\, gigiH 1 ,with i = 1, 2, 4 are the sum of all the trees with only 
end-points at scale < and with propagators g^\ see (2.11), and in which the factors Z g~ 1 , k > h 
are replaced by 1 + z\. 

3) The terms contributing to Phi&h are by definition independent from g\. kl k > h. 

4) The terms contributing to G l hl G 2 h , G\ by definitions depend linearly from gi_ k , that is they are 
vanishing if g ik = for any k and their second derivatives respect to g\, k are also vanishing, while 
the first derivative are not vanishing. 

5) The terms at least quadratic in gi are included in J2 k k , gi,k9i,k' k k , and by (3.42) it holds 

\Hl Kk ,\<Cv hl e ^^ h -^ (4.9) 
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6)In R£' we include; terms depending from Vh or 8h\ terms with at least a propagator r^(x — y), 
see (2.11); or terms with at least an endpoint at scale 1. 

Note that the above decomposition is obtained by an analogous decomposition over trees, so that 
the determinant bounds of §3 are still valid. 

In writing (4.8) we have used that the beta function contributing to gi has at least a gi, in fact 
consider a contribution to the antiparallel part of gi; it is not invariant under the transformation 
— ► e ±a ijf^ a and ip^ 1 a — ► ip^ 1 a while the terms corresponding to gi and 54 are invariant. 

The flow given by (4.8) is very difficult to study; luckily dramatic cancellations appear, given by, 
if g h = max fe > h (|^| + \g\\ + \g%\) and p, h = max fe > h (|^| + \gl\), the following result. 

Theorem 4 (Partial vanishing of the Beta function) The functions 0^, (i\,G\,G A h ,G\, 



for \vh\ < s are such that, for a suitable constant C 

\p 2 h (»h,..,» h )\<Cj2fa eh \/3 4 h (» h ,..,» h )\<Cjifa eh (4.10) 

\Gl(g h ,..,g h )\<Cgfa eh \G£(g h , ..,g h )\ < Cgfa eh (4.11) 

\Gl(g h ,..,g h )\<Cgfa eh (4.12) 



The above lemma says that a dramatic cancellation happens in the series for the above functions; 
each order is sum of many terms 0(1), but at the end the final sum is 0(-y eh ), that is asymptotically 
vanishing. We call such property partial vanishing of the Beta function (partial becouse the O(gf) 
terms are not vanishing). 

Assuming the above lemma, which will proved in the following two sections as consequence of 
suitable Ward identities, we can prove that the flow is bounded for any gi t o > 0. 

def def 

We proceed in the following way. We first assign a sequence v_ = {vn}h<i, S= {Sh}h<i not nec- 
essarily solving the flow equation for v, 6, but such that |^|,|<5h| < cll-f eh , f° r an y h < 1. We 
then solve the flow equation for gi : h, parametrically in v, 6, and show that, for any sequence v_,8 
with the supposed property, the solution g{v_,8) = {ffVi(l£, <!), .92,h(l£, <£), g4,h{Ki <L)}h<i exists and 
has good decaying properties. We finally fix the sequence v_ via a convergent iterative procedure. 

Lemma 1. Assume that \v h \, \S h \ < cUj eh f or an U h. For U > and small enough the flow is 
given by, for any h 

\92.h -.92,0 -.9i.o/2| < C/ 3/2 \94,h - 94.o\ < C/ 3/2 0<ffi, fc <- ^ (4.13) 

1 — a/ogi.on 



Proof. By using that \i/h\, \Sh\ < cll-f eh it holds that 

\R\\ < CU 2 -f 0h (4.14) 
It is convenient to introduce #2, h — %92,h — 9i,hl then using (4.10) and (4.14) 

§2fi-i=h,h + ^D htk + ^(2I% tk -Dl k ) + ^g ltk g ltk ,H hiktk ,+R h (4.15) 

k>h k>h k,k' 

with 

Dh.k = Pl(Vh, ■••M/nMfc,Mfc+ii -,Mo) - PliPh, ...Hh,Vh,Vk+i, -jMo) (4-16) 
D h,k = G h(9h, ■■■9h,g' k ,9k'+i, -,go) - Gl(g h , ...g h ,g h ,g k , +1 , ..,g ) i = 1,2 
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and a similar equation for g^h', Hh,k,k> verifies (4.9), Rh (4.14) and 

\D h , k \ < C^-Vu^h-Q^ \Dl k \ < CU^-^gu-g^ (4.17) 

Assume that for k > h 

0<(H,k-i<z " 1 \ , y. \9k-.-9k\<[uh ek + \ g ).° J 2 ] (4.18) 

1 - a/3gifl(k - 1) I - a/3g lfi k 

We have then to prove that such inequalities hold for k = h — 1. Noting that 



k=h k—h k—h 



we obtain 



Moreover 



J29i,k9i,k>H 2 hXk ,<CUgl h (4.20) 
fe.fe' 



E < E cc^ 7 -^*- fc ) E i^V* + f ^;^. ] 2 ! < 

fe>h fe=-l fc'=/i ' 

C 2 C[/([/t 7 ^+ V |fc~% 2fl ^- fc )[ i g ^ ° I 2 ) (4.21) 
k ~_ l I - a/ igifik 

and the last addend can be bounded by 



Then by (4.15) we get 



\h,h-i - h,h\ < C 3 (U 2 7 eh + U( 9l / ) 2 ) (4.23) 



and 

k—h k ' 

In the same way in the flow for 54 we use that there are no second order contributions quadratic 
in gi h- Finally we write, using (4.8) and the short memory property (namely that ^ e ( h ~ k )g 1 k < 
CgiM) 

9\,h-\ ~ 9i,h < ~-^9i,h9i,h-i (4.25) 

or 

gi.H-1 < -rfr— (4.26) 

1 + 3^1, h 



and as is an increasing function and by induction < g\_h < 1 _l 1 g ° ah so that 

gi,o(i- ffegLo)" 1 91.0 27) 

yhh ~ 1 - 1 + 1^0(1-1^0)-! - l _a gi(j{h _ l) ■ ' 

Moreover gi^-i = ffi./^l + 0(£/)) by (4.8), and g 1Jt > so that gi. h -i > 0. ■ 
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4.3 The choice of the counterterms. 

In the previous section we solve the flow equation for parametrically in any sequence v_ = 
{vh}h<i, §_ = {Sh}h<i such that \i>h\ < cU-f eh , \o~h\ < cll-f eh f° r an y h. We show now that indeed 
we can choose v, 5 so that v_ = {vh}h<i, 5 = {<>h}h<i verify such a property. 

Lemma 2 There exist sequences v_ = {vh}h<i, H = {5h}h<i such that \i/h\ < cUj eh , \Sh\ < cll-f 0h - 



Proof It holds that 



/f^C+C ( 4 - 28 ) 



where us given by a sum of trees with no end-points v k , 5k and only propagators g\ (2.11); 

by the symmetry in the exchange x,xq of g^, and remembering that (3^ — J2 t [z(t) — a(r)] it 
holds that 

\Ps k !\ < C'U 7 26h (4.29) 
A similar decomposition can be done also for 

^=^1+^ (4-30) 

again with 

„28h 



ti\<CU-i Mh (4.31) 

by the parity property g\ w (x, y) = —g\ w (y, x). If we want to fix v, 5 in such a way that u_ oa = 
5-oo = 0, we must have, if (fi,5i) = {v, 5): 



l 



v = - ^ 2 P k v{gk,5 kl v k ;...;g ll 8 l ,v l ). (4.32) 



k— — oo 
1 



S=- Ps(9k,$k,vk; ■ ■ ■ ;5i,<5i,fi) ■ (4-33) 



k— — oo 



Note that in (4.32), (4.33) g k = g k {v_A)- 
If we manage to fix v, 5 as in (4.32), (4.33) we also get: 



»h = ~ E 7 fc " h " 1 ^(flfc, 8 k , vk\ ■ ■ ■ ; <?i, • (4.34) 



^ = -^^(fffc,^,^;---; gi,&i,vi) ■ (4.35) 

Let 9Jlg be the space of sequences v = {f-oo, • ■ • , vi}, 5 = {5-oo, •, ^i} with small || • ||g norm, 
namely the space of sequences v_ satisfying: 

141 < Y k ,Wk\ < i ek 

We look for a fixed point of the operator T : Tie ~ f 9Jle defined as: 



T(yh) = -J2 l^-^Kdk&v), v k -...- g^v x ) . (4.36) 

k<h 

T{5 h ) = -Y J Ps (9k {5,v),5 k ,v k ;...;g 1 ,5 u v 1 ). (4.37) 



k<h 
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First note that, if U is sufficiently small, then T leaves Tig invariant: in fact 

|(Tv)/,| < ^2 2c 1 U-i 0k -f k ~ h < cU-f eh \(TS) h \ < ^2ciC/ 7 efc < cU~f 0h (4.38) 

k<h k<h 

Furthermore we find that T is a contraction on Tig: in fact 

\(TS) h - (TS') h \ < ]T \Pli9kia, S),u k ,S k ;...)- fi{g k {j/,#), v' k ,^...)\< 

k<h (4.39) 

<c"U^ h {\\v-!/\\e + \\S-S!\\e} . 

and a similar equation holds for v. Then, a unique fixed point v*,5 for T exists on Tig. ■ 

By the above Lemma we have found S(t,pF, U), v(i,pF, U); inserting them in (2.1) and using the 
implicit function theorem we get Pf{U, /z), t(U, ji). 
Finally from an explicit second order computation we obtain that 

Zh = a[gl h + gl h + g\ h \ + ^ (4.40) 

with a > is a suitable constant, and using the previous results on the flow of g^h, Vh,5h we get 
lim^oo ^ = 1, where r, = a[gl_ x + gt_ x ] + 0{U 3 ). 

5. The reference model and proof of Theorem 4 



5.1 The model. 

In order to prove the partial vanishing of the Hubbard model Beta function expressed by (4.10), 
(4. 11), (4. 12) we introduce a reference model written directly in terms of grassmann variables, with 
an ultraviolet cutoff and an infrared cutoff j h with linear dispersion relation and in the continuum. 
We study the reference model by Renormalization Group and we show that the Beta function of 
this model is asymptotically vanishing as a consequence of Ward identities due to the formal local 
chiral gauge invariance (which is however broken by the presence of cutoffs); then we prove that 
the Beta function of the reference model coincides partly with the Beta function of the Hubbard 
model, so that we can deduce the partial vanishing of the Hubbard model Beta function from the 
vanishing of the reference model Beta function. 
The partition function of the reference model is 



/ 



P L (dip)e VL (5.1) 
where the propagator is 

&,^-y) = ji E <sJw5sS (5 - 2) 



with CrJ(k) = EL-oo fHh, k) and 



,/3/2 r-L/2 

Vl = 12 / dX ° / dx [92^t,u,a^,u,At-u J ^.-u.a + 

,. J-0/2 J-L/2 



1-0/2 rL/2 

dx 

-13/2 J-L/2 

g2^tv,<j^,u,a^i-u-a^-u-a + ^x^^x^^x," - CT ^x, w -a] ( 5 - 3 ) 

Note that the model is not SU(2) invariant, as the interaction depends from the spin if g% ^ g\. 
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The Grassmann integration can be done by a multiscalc analysis essentially identical to the one 
described in §3; however the symmetries of the interaction imply that the local part of the effective 
potential (3.29) is replaced by 

/-/J/2 rL/2 

CV i = 22 dX ° / ^[^jV^aC^^-^aC,^ 

V J -W J-L/2 

(5.4) 

+9i,jbPt^,^2^,^t^-^^-^ 
Note in fact that the analogue of Vh,5h are vanishing by (in the limit L,[3 — ► oo) parity and 
invariance in the exchange (x,x ) — > (xq,x) ; moreover the reference model is invariant under the 
total gauge transformation ^ u a — > e ±a "-°i/'J u CT for any values of a W)0 -, so that terms of the form 
' l Pi,(T' l P~-Lu,cr' l f J -Lj,-o' l l , Z.-(j cannot be generated in the integration procedure as they violate such 
symmetry. Note also that, due to the compact support of the cutoff in (5.2), the running coupling 
constants at scale k > h of the theory with infrared cutoff 7^ or are identical. 
It is easy to verify that a tree expansion similar to the one described in §3.5 holds also for the 
reference model, and that the analogue of Theorem 3 holds also in this case. We will prove in §6 
the following Lemma. 

Lemma ZAssume that g — maxQg^l, Isfl, Ifi^l) * s sma ^ enough; then for any integer j < 0, in the 
limit h — > — 00 for a suitable constant C 

\~9h-9°2\<Cg 2 \%j-<$\<C? \~g4 tj -9 4 p \<Cg 2 (5.5) 

Moreover 92,ji92,ji9i,j have a limit as h — > —00. 

It is an immediate corollary of Lemma 3 and Theorem 3 that = {g^ fc> 9% fc> 9i,k) are analytic 
functions of v\ = (g^, g\, 94) around (0,0,0). Note that analyticity in the coupling around the 
origin holds for the reference model and not for the Hubbard model. Finally the analogue of the 
flow equations (4.8) is given by, if v% = (ff^fe) 92 fc> 9i,k) 

-glj-^gh + ^ivf,:,^) (5.6) 

94J-1 = 9i,j + -,Vq) 

We can rewrite the above equations as, for j > h 

9h-i-9h+^°(vf,..,vf)+J2D^ 

k>j 

k>j 

~94j-i = 94J + vf) + ]T Df° 

k>j 

with,for a = {2,6), (2,p),4 

D "k = -A) ~ -^f,vf,vf +1 , -,4) (5-8) 



5.2 Vanishing of the reference model beta function. 

The Beta function is an analytic function of vf and it can be written as, if a — (o, 2), {p, 2), 4 

Pf{vf,...,vf)= x: ^^^(^r^-)" 2 ^)" 3 
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We define n = n\ + n 2 + n 3 and n = (m, n 2 , n 3 ). Note that 



6" =6" +OfT e - J ) 

"j,ni,ra2,ra3 11,12,13 ^ I / 

Consider ^j jnii „ 2i „ 3 and 6^ n with k < j; for any tree r contributing to /3fc there is a tree 
contributing to (3j; in fact we can perform a change of variables in the propagator g l (k) respectively 
k — > 7 J k and k — > 7 fc k, so that in one case the propagator is /(7~ I+fc k)Z)J 1 (k) and in the other 
/(7~ l+ -'k)-D~ 1 (k) = /(7~* +fe 7 _fe+ - 5 k).Dj 1 (k); hence for each tree contributing to bj there is a tree 
contributing to bk, in which the scale of each vertex is shifted by j — k; there are extra contributions 
to b k with at least a vertex with scale > k — j; such trees have the root at scale k so that, by the 
short memory property, bj — bk = 0(7 with < 6 < 1 a suitable constant, and taking the limit 
k -> -oo we get bj = b + 0(-y ej )- 
We will prove the following result. 



Lemma 4 Assume that Lemma 3 holds; then for any (ni,n 2 ,n 3 ) 

K u n 2 ,n 3 = (5-9) 



Proof. The proof is by contradiction; assume that, for some n = (rti,n 2 , n 3 ) with n — ni+n 2 +n 3 

f3j(vt, ...,#) = Knighri^r^r + o<$?+\ (5.10) 

with b= a non vanishing vector, and that for all ri\ + n 2 + n 3 = n < n — 1, bft is vanishing. From 
Theorem 3 and Lemma 3 vf are analytic functions of v\ — (g%, .g 2 , .94), that is 

«t = «1 + E 4 j) (52) ni (^)" 2 (54)" 3 + 0((^) S+1 ) (5.11) 

n<n 

and for any fixed j the sequence cl is a bounded sequence. Inserting (5.11) in the Beta function, 
using analyticity and equating the coefficient of (g 2 ) ni (g% )™ 2 (ff4)" 3 with m + n 2 + n 3 = n < n — 1 
we get 

4 j_1) =4 j) + E ^ + 0(7^) (5.12) 
fe=j+l 

where the last sum represents the contribution of Dj^, so that 

\4 k \< 7 ~ e ^D n sup (5.13) 

2<m<ri-l 

where we have used that Djj~ is at least quadratic in the running coupling constants, and D n is a 
suitable constant (in j). Note that 

\4~ l) ~cn\<Dn E t E SU P l^-^l+l^")] (5-14) 

., 7-/11 2<m<n— 1 

The above inequality implies by induction that, for n < n — 1 

sup |cf - ck| < C" 7 f * (5-15) 

2<m<n 

for a suitable C; assume in fact that it is true for & > j and from (5.14) we get 

- ck| < C^Dn IE 7- e(fe - / H7 |fe +7 fj '')+7 ej ']<^C' n - 1 5 S 7 f0 '- 1) (5.16) 

j' = -oo fe=j' + l 
2//c66raio/2008; 13:13 23 



so that (5.15) holds for j — 1 if C > KDn. On the other hand (5.15) implies 

< C n ji W" 1 ) (5.17) 
Writing now the analogous of (5.12) for n = n we get 

c i" 1) = c i ) +^«+^ (5.i8) 

which can be rewritten as 

c t~ 1) = c i' ) +^+°^ fj ) ( 5 - 19 ) 

so that crJ is necessarily a diverging as j — > oo, and this is a contradiction. ■ 

5.3 Partial vanishing of the Hubbard model Beta function (Proof of Theorem 4)- 

We compare the Beta functions of the reference model with the function appearing in the flow 
equations of the Hubbard model. 

A) Let us start considering first the reference model in the spin symmetric case, that is if <?2.o = 9% o- 
In such a case by the same arguments used in Appendix , for any k g% k = g v 2 k and 0% k = 0^ kl so 
that the flow equation (5.6) reduces to 

92,h-\ = <?2Ji-l + Ph{g2,h,9i,h\ ",52,0,54,0) 

54,/i-i = .94J1-1 + Ph(92,h,g4,h', ■•; 52,0, 54,0) 

It holds that the functions 0\ and 0\ essentially coincide with the functions 0f, 0^ of the Hubbard 
model defined in (4.8); that is, if fih — (g2,h, 94, h), for a suitable constant C 

0Z( f i h ,..,» h )-f3t(» h ,..v h )\<Cf4 ll eh (5.20) 

\^{fjL h ,..,H h )-^{fjL h ,..^ h )\ <Cnfa eh (5.21) 

The above equations prove (4.10). In order to prove (5.20) and (5.21) we note that by definition 
the only difference between 0^,0^ and 0^,0^ is that in one case the model is defined on the 
continuum and in the other case on the lattice. In momentum representation this means that the 
delta functions in (3 are defined as L06kfi^k o ,o while in [3 are defined as in (2.16). The difference of 
the two delta functions slightly affects the non local terms on any scale, hence it affects the beta 
function; however, it is easy to show that this is a negligible phenomenon. Let us consider in fact 
a particular tree r and a vertex v £ t of scale h v with 2n external fields of space momenta k' r , 
r = 1, . . . , 2n; the conservation of momentum implies that Y] " 1 o~ r k' r = 2irm, with m = in the 
continuous model, but m arbitrary integer for the lattice model. On the other hand, k' r is of order 
j hv for any r, hence m can be different from only if n is of order j~ h " . Since the number of 
endpoints following a vertex with 2n external fields is greater or equal to n — 1 and there is a small 
factor (of order fi h ) associated with each endpoint, we get an improvement, in the bound of the 
terms with \m\ > 0, with respect to the others, of a factor exp(— C-f~ hv ). Hence it is easy to show 
that the difference between the two beta functions is of order n\^ eh . 

B) In order to prove (4.11) we consider the reference model with g®^ ^ g% , so that there are 
three independent running coupling constants. We have seen that, for a = (2,p), (2, o), 4 

ft{vf;,..vk)= £ ^^^.^feM-MffL]" 2 ^,,]" 3 (5.22) 

ni,n 2 .n 3 

On the other hand we can write the functions G% (4.8) in the Hubbard model, a — (2o), (2p), 4, as 

Gt = E <t,i,n 2 ,nMA[92,h] m2 M m3 (5-23) 

m2,m 3 
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The coefficients 1 n2 n are given by sum of trees (or product of trees, for the presence of the 
z\ terms) with (in total) one end-point g\, mi end-points gi and mz end-points g^\ the SU(2) 
invariance of the Hubbard model implies that G|° = G 2p . To g\ and gi correspond two terms, 
the parallel or antiparallcl part, see (3.35), and we can associate to the endpoints of the trees 
contributing to 1 m2 ms an extra index distinguishing the parallel or antiparallel part; then we 
can write 

C h,l,m 2 ,m 3 = C h,m°,m^,m°,ml,m 3 (5-24) 

It holds that 

C h,l,m 2 ,m 3 — E C h,0,l,m°,ml,m 3 (5.25) 
m^-\-m^ — m 2 

that is only the spin parallel part of gi can contribute to G\ or G\; in fact making the the global 
gauge transformation ipf^ — > e lcr ^ (J and ip^ a — > V^i o-i the antiparallel part is not invariant, 
while the spin parallel (and the gi, 54 interactions) are invariant. 

Finally note that the spin parallel gi interaction is equal (up to a sign) to the spin parallel gi 
interaction, so that, for a. = (26), (2p), 4 

C 0,l,m°,m£,m 3 = "C>^Lm 3 = (5.26) 



C)It remains to consider (4.12); we can consider cquivalcntly the contribution to the spin parallel 
or the spin antiparallel, as they are equal by SU(2) invariance of the Hubbard model, that is 
G l° = Gl h- We consider the spin parallel part and we can write 

m2,m 3 

with 

r lp = V V r lp 

L -l,m 2 ,m 3 m°,m.P,m°,mP,m 3 

The single gi interaction cannot be antiparallel, again because making the global gauge transfor- 
mation ipf^ — » e la i/)f a an d $-10 ~> Tp-i ct' tne antiparallel part is not invariant, while the spin 
parallel (and the gi, c/4 interactions) are invariant. Hence 

C l P m 2 ,m 3 = E C ? l,m 2 °,m^m 3 (5.28) 

and 

u 0,l, m^,m^,m 3 u m'^,m 1 ^ + l,'in 3 

as the contribution (Ip) and (2p) are identical. 



c lp P =b 2p p., =0 (5.29) 



6. Ward identities for the reference model: Proof of Lemma 3 



6.1 Dyson equations 

Let us now prove (5.5), extending the analysis in [BM2], [BM3], [BM4] to the spinning case. We 
derive a number of Dyson equations relating some Schwinger functions of the reference model. Let 
us start from, if p PiWi<T = jfc £ k Vk.^k-pA^ 

< V'k 1 ,+, CT V'k 2 , + , CT V'k3,_,_ CT V'k 4 ,_,_ <7 >t= 

3 -(k 4 ){G^(k 3 )[ 32 < V'k 1 , +;a V'k 2> +, CT Pk 1 -k 2> +,a >T 
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+92 < V'k 1 ,+, CT V'k 2 ,+, CT Pk 1 -k 2 ,+,-a >T 
+94 < i> + ll+ ^k 2 , + ,aP^-^2,-,-^ >T} 

+ J dp[g% < Vk li+;a Vk 2 ,+ >CT V'k3 > - >CT V'k 4 -p > - > - <T Pp,+ > a >t + (6.1) 

J dpgP < Vk 1 ,+, CT V'k 2 ,+,aV'k3,-, CT V'k 4 -p,-,- CT Pp,+,-T >T + 
y rfP54 < V'k 1 ,+,aV'k 2 ,+ > aV'k3 > +, CT V'k 4 -p,+,- CT Pp,-,-a >T]} 

where 

^(k) =< >^ (6-2) 

Similar Dyson equations holds for < ^£ >+iff V£ 2 , + ,^,_, ff ^,-,a >t and < ^£,+,^,+,^3,+,^, 
The Renormalization Group analysis of the preceding sections easily implies (for details, see [BM3]) 
that, if |kj| =7 fc ,t= 1,2,3,4 

< C+^+X,-,-^,-,-. >t= Gt >( ,(ki,k2,k3,k4) = 7- ih Z^[g° 2 , h + O(gl)} (6.3) 

if g~h = su Pfc>/idS2 fel + Iff fcl + 1 54 1 fcl)- ^ n t ne Dyson equations appear the functions 

G t]a,u>', a ' ( k - q, k, q) =< V'k.u-^V'q^aPk-qX^' >T (6.4) 

G + 1 ^' CT '( k i' k 2>k3,k 4 -p;p) =< i> + u+}a ^ 2t+ ^ + 3 _ } _ a ^ 4 _ P} _ t _ a P P ,u,>,*> >t (6.5) 

Either such functions or the Schwingcr functions can be obtained by deriving the Generating 
functional 

w(<f>, j) = log J p i (#) e - ViW+ ^. . / ^( J -.-.^i.«.^-.u...+<u-.^x. w ,.+< w ,^, w .j (66) 

with respect to the external fields J x , w ,<t or (p^u a- 

The functions G 2,1 ,G 4,1 are related by remarkable H^ord Identities to the Schwinger functions 
G 2 ,G 4 . In fact, by operating in (6.6) the (local) Gauge transformation ip+ a — > e ±ax xp+ a , 

ip+.-a — » V^-a and V^iff — » V'-.itr and deriving with respect to cj>+ +M ,cj>+ +IT , we get, pass- 
ing to momentum space, the following Ward Identity 

D + (p)G^ 1 a , +;a (k 1 -k 2 ,k!,k 2 ) = G 2 + (k0 -G 2 + (k 2 ) + A^ i+;CT (k! -k 2 ,ki,k 2 ) (6.7) 

where 

A !/^>' ( k - q, k, q) =< V'k.u^V'q.^^Pk-qX,^ >T (6.8) 

8pk-w = j dk 'G(k', k' - q^^C'-k+q,^ (6-9) 

G(k+,k-) = (C fci0 (k-) - l)D w (k-) - (C h , (k+) - l)D w (k+) (6.10) 

The cutoff function Ch,o in Pi, (dtp) destroys the local Gauge invariance of the theory, and it is re- 
sponsible of the correction term A 2 ' 1 in (4.3). As explained in §4 of [BM2], AW) = G(k, q)#W (k)pW (q) 
is non vanishing if at least one among i or j is or h; this means that either at least one field in 
Sp is contracted at scale 0, or at least one field in Sp is contracted at scale h. We can split the 
correction term in the following way 

A 2 ' 1 , , ^A 2 ' 1 '", , + A 2 ' 1 ^, , (6.11) 
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where in Aj ^™ , CT , there are all the contributions with one of the fields in dp contracted at scale 
0, and A 2 '^, a , is the rest. It is easy to check that 

|A^t>'(ki -k2,ki,k 2 )| <9h^- (6.12) 

This follows from the bound \ \ < 7 fe ~ J 7 z and noting that the factor gives the correct 
power counting for the marginal terms linear in J, see [BM2]; note also that the contributions of 
order in the cancels out. 

The analysis of A 2 ' 1 ^, a , is more complex; there are other remarkable identities (first discovered 
in [BM3] for the spinless case) called correction identities to the functions G 21 . It holds in fact 
the following Lemma. 

Lemma 5 There exists functions Vv,±a such that \v Ut ± a \ < Cejh and 

A^ +;CT = < iCTJ D + (p)G 2 ^ i+;CT 
+^_ rT D + (p)G 2 ^_ rT + < CTJ D_(p)G^,_ i(7 + iS^D-WGQ^ + (6.13) 



with, if |ki| = |k2 1 = 7 ,! 



\H 2 ^%^ - k 2 ,ki,k 2 )| < (6.14) 

h 



for some constants C and < 8 < 1 . 

The above identity says that the correction A 2 ' 1 '" can be written in terms of the functions G 21 , 
up to a term which is smaller than 0(Y' 1 )- We will call H 2 ' 1 = H+^ a + r7 + , and 



\B^\Vi-V 2 MM)\<g h \-. (6.15) 

By the phase transformation ip^ _ a — » e ±a *V+ _ cr , CT — * ^+ <j an d ±<r ~~ * ± CT ' an< ^ usm g 
a correction identity similar to (6.13) we find 

-< a D + (p)G^ CT;+iCT +(l-<_JZMp)^ 

(6.16) 

where H 2 ' 1,a verifies a bound similar to (6.15). 

In the same way by the Gauge transformation ip^; a — > e ±ax V^ ^ _ CT — » ^ _ CT and ±cr — ► 
^+ ±cr > an( i usm g a correction identity similar to (6.13) we find 

(6.17) 

where H 2 ' 1 ' 11 verifies a bound similar to (6.15). 

Finally by the Gauge transformation f/' d :,_ (T — > e ±Q!x V- -cr> V- cr — * V 1 - u an d ^+ ±cr — * ^+ ±<r we 
get the following WI 

-^D^p^^-^D^G^ 

(6-18) 

Bounds like (6.15) and (6.14) hold also for H b ' ,H 2 ' 1 ,H d ' . It is easy to see from some algebra 
that the above relations imply 

D + (p)G 2 ^ a (k a - k 2 , ki, k 2 ) = G^kx) - G+(k 2 ) + (1 + FDHl' 1 + FfrH 2 ' 1 + F^H 2 - 1 + F}H 2 / 

(6.19) 

D^Gl^Jk, k 2 , k 1; k 2 ) = (1 + F 2 )/^ + ifff 2 ' 1 + F 2 !! 2 ' 1 + F 2 H 2 / (6.20) 
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D_ (p)G 2 ^_^i k 2 , k u k 2 ) = (1 + F^)H 2 ' a + FfH*' 1 + F?H^ + FjH 2 / (6.21) 

ZMp)G",-,-.(ki - k 2 , k l5 k 2 ) = (1 + Ft)Bl' a + F b 4 fff + F^' 1 + FjH 2 / (6.22) 

with F a ,Fi, 7 F C7 Fd are combinations of the v, with the property that if \uf \ < Cg, then \Fi \ < Cg~h- 
Then (6.19) really provides a relation between G 2 ' 1 and G 2 up to bounded corrections. 

6.2 Proof of Lemma 5. 

We introduce the generating function for H^' 1 

J P L (dij)e- VL+Tl+ f dkdpJp "+-« D+{p) ^t+.^ k +^+° (6.23) 

e - f dkd P j p [^_ a i? + (p)v>+ +i _^k +p , + ,^ 
where 

T\ = y dkdp J p G(k, k + pX + , CT Vk +P;+ , CT (6.24) 

The analysis proceeds essentially identical to the one of §4 of [BM2]. After integrating the ip° 
field, we get in the effective potential a sum of monomials of the form WJ m ipi...ip n ] we extend 
the definition of C to monomials of this kind by requiring that it acts non trivially only on the 
terms linear in J and quadratic in ip, as a power counting argument shows that they are the only 
marginal terms. 

Consider now the terms in which Ti is contracted; they are of the form 

E/<*p/ dk+J p^^ k+ -p,^^^^^^ (6-25) 

LJ.CT 

where F^ ~+\ & 5 i s given by all the terms obtained contracting both the ip fields in Ti while i*i 



are is given by the terms obtained leaving external one of the tp-Helds of T\ Both contributions to 
the r.h.s. of (3.3! 
it is of the form 



the r.h.s. of (3.39) are dimensionally marginal; however, the renormalization of f[J\ is trivial, as 



F^i (k+ k") = [ I^oOO - a^gp^gg) - "°( k+ ) G P)( k +) (6.26) 

D + [k k ) 

or the similar one, obtained exchanging k + with k . 

By the oddness of the propagator in the momentum, (0) = 0, hence we can regularize such 
term without introducing any local term, by simply rewriting it as 

<- + >+, k -) = = Zt+ ] -lr " Mk+) [G(2)(k+) - G(2)(0)] • (6 - 27) 

As shown in [BM2], by using the symmetry property 

ff«(k) = -^)(k*) , k=(fc,fc ), k* = (-fc ,fc), (6.28) 
' ^ - ~ can be written as 

F 2^.*.cA k+ , k ~) = -TTT^ Mo,^,*(k+,k-) +piAi ia)i<ri0iff (k + ,k-)] , (6.29) 
where A ijUJ ,<7,u,a(k + , k~) are functions such that, if we define 

£F^l,°fi,±° = 7rrTbo^o,+, CT> <;,±a(o,o)+p 1 A 1;+iCT;(;i±CT (o,o)] , (6.30) 

L> +\P) 
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then, 

£^,^, ±<T =^(p)^!f' ±CT , (6.31) 

where z 3 :^' ±(T are four suitable real constants. 
Consider now the terms in which the v u ^ a arc contracted; we define the localization operator on 
such terms as 

C J dkdpD u (p)W-iAK^-p)^ ±!7 ^ +P ^±a = J ^pA,(p)^L(0,0)VC,±>k +P ^±. 

( ° 2) 

We define v-\ yU j,±a = Zjf' a + W Wj±cr we get that the local terms linear in J arc 

J dkdpJ p [^ 1 ^ a D + (p)yj+ + ^ +p+a + (6.33) 

i£i, +> _,I> + (p)^, +i _,C^ +> -^ 

We can iterate the above procedure; at the integration of the generic scale the terms quadratic and 
linear in J in the effective potential are obtained contracted a Ti vertex (in such a case one of the 
two fields of T\ is necessarily contracted at scale 0) or a Vk,u>,a vertex; in both case the preceding 
analysis can be repeated and the local terms linear in J are, for k > h 

J dkdpJ P K ;+ ,^+(p)< + ,^k +P , +ff + (6-34) 
<+,_ ff I> + (p)Ch_,VW 1 _ < , + ^ 

We have then obtained an expansion for H^ 1 in which new running coupling constants ap- 
pear, namely Vk.u.±<r', the analogue of Theorem 3 ensures convergence Vk, u ,± a are small or any 
k > h. The beta function for Vk,u,±a has the following form 

Vk-\,u,±o = V k , + ,±a + /3i;± CT (Wfc , -,Vq) + Pl% a (Vk, Vk-,VQ, v ) (6.35) 

where by definition 0^±„{v^, Vk-—, v o , ^o) is obtained contracting a Vj while ^'^(v^ , Vq ) is 
obtained contracting Ti and 

\fc k ±M,-A)\<c- 9 ki ek 

for some constant < 9 < 1 . The presence of the factor -f ek ln tne above bound is due to the fact 
that, for the support properties of the function C(k + ,k~) discussed after (6.24), one of the fields 
of T\ is necessarily contracted at scale 0. 

In fact we can show (proceeding as in the proof of Lemma 3, or in §4.6 of BM3) that there exists 
a sequence v k such that \ fk,ui,±a 

| < CghY k by solving 
k 

vk,»,±o = - ]T -A) + ^o)} (6.36) 

k'=h+l 

This shows that there exist Vu,±o such that Vk,u>,±a — 0(j ek ). 

We have then find an expansion for H^iki — k2,ki,k 2 ) very similar to the one of G 2,1 , but 
in which each tree contributing to iJ 2;1 (ki — k2,ki,k2) have an extra r ) 6h ] in fact or there is an 
endpoint v k (and we use that Vk,uj,±a = 0(-f 6k ) and the fact that, as the dimension are negative, the 
value of the tree has an extra ^ Bi - h ^ k )) or there is an endpoint Ti contracted at scale (hence, as 
the dimensions are negative,the value of the tree has an extra 7 6 " 1 ). ■ 

Inserting (6.20) in the Dyson equation, and using (6.15) and (6.14), we see that the first three 

— Ah 

addends of the Dyson equation are given by (g® + 0(cj^))' 1 -^- . 
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We have to consider now the last three addenda in the Dyson equation; let us start by 

/ dp[g° 2 < ^u+JL+^L-.-^-p.-.-aPp^ i • :,l - :i7: ' 



Let us call 

-»4,1 



G+ !(TA/(T ,(ki,k 2 ,k 3 ,k 4 -p;p) =<i>i u+ ^^ + ^i 3 ,-,-^ i - v ,-,-aPv^'^ >t (6.38) 

As |k 4 | = -f h the support properties of the propagators imply that |p| < 7 + -f h < 2j, hence we 
can freely multiply G+ 1 in the r.h.s. of (6.37) by the compact support function xn{l~^ m IpDi with 
jm = [1 + l°g 7 2] + 1. It follows that (6.37) can be written as 

J rfpXAf(p)G^ 1 (p;ki,k 2 ,k 3 ,k 4 - p) + J dpxM(p)G+ 1 (p; ki, k 2 , k 3 , k 4 - p) (6.39) 

where xm(p) is a compact support function vanishing for |p| > ^ h +i^ 1 and 

Xm(p) = £ AM- (6-40) 

h p =h+j m 

Note that the decomposition of the p sum is done so that Xm(p) = if |p| < 2j h . It is easy 



to show that the first term in (6.39) is bounded by 0{^ b -L i — ), see [BM3]. Regarding the second 



h 



addend we will use the following Ward identities 
(l- V a + . a )D + {p)G? M V %,_ a D + ( P )G^ + ,_ - < CTJ D_ (p)G^_, a v\_„D_ (p)G« = 
G* (kx - p, k 2 , k 3 , k 4 - p) - G% (k 1} k 2 + p, k 3 , k 4 - p) + i#i( kl> k2j ks) k4 _ p; p) ( 6 .4i) 



-^D^G^^+a-^.jzMp)^^^ 

(6.42) 

-^ a D + (p)G^ + , a -^_ a D + (p)G^ + ^ 

(6.43) 

-<,IMp)^, + .a-<-«,0+(p)^^ 

G 4 (ki, k 2 , k 3 - p , k 4 - p) - G 4 (k l7 k 2 , k 3 , k 4 ) + ^''(ki, k 2 , k 3 , k 4 - p; p) (6.44) 

where the functions Hf' 1 are defined in an analogous way to the functions Hf' 1 . It is easy to see 
from some algebra that the above relations imply 

Z? + (p)G^_ CT (k!, k 2 , k 3 , k 4 -p; p) = (1+G?)[G 4 (ki-p, k 2 , k 3 , k-p)-G* (ki-p, k 2 , k 3 , k-p)] 

+G a 2 [G% (ki, k 2 , k 3 - p, k 4 - p) - G 4 (ki, k 2 , k 3 , k 4 )] + (1 + Gl)H A a ' a + Gtrf' 1 + G 5 a H^ + G 6 a H 4 / 

(6.45) 

with G l a = 0(g). From (6.19) we can decompose / c?PXm(p)G 4 J 1 as sum of several terms; the one 

involving G 4 (ki, k 2 , k 3 , k 4 ) is vanishing while the other three terms involving the other functions 

-2-3/1 

G 4 have a bound 0{^jp. — ), see [BM3]. Finally the following results holds 
Lemma 6 // the functions Vu.± a are the same as in Lemma 7, it holds that, for i = a,b,c,d 

if 4 ' 1 92 ~ 3/l 

dp9AkA) Dt(p) l - C Zf~ (6 ' 46) 



4.1 
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Inserting all the above bounds in the Dyson equation (4.1) computed at momenta |k,| = j h ,i = 
1, 2, 3, 4 we have completed the proof of Lemma 3. 



Lemma 6 is proved considering 

Gt(ki,k 2 ,k3,k4) = — - — — + — - W 



(6.47) 



where 

W = log J p(d^) e - T (^)+ ,y i aT i(^)+ ,y 2r2(^)+^ a T 3 (^)+^T4(^)+X; ij / dxttl^-^+t+^-J i (g 48 ) 

T W0 = Tj3 E * M{p) If3 E C+ ^p) ^ (V i k, + ,^k-p, + ,.)^ 4 - P ,-,-.^g-(k4) , (6-49) 

P K 

Tl M = ^ExM(p)^E^k, + >k- P , + >k 4 -p,-,-^k 4 3-(k 4 ) , (6.50) 

W) = ^E^(P)^E^k )+ ,-^k-p, +> -JV i k 4 -p > - > -^k 4 5-(k4) , (6.51) 

w) = ^E^(p)^E^^^k,-,^k- P ,->k;-p,-,- CT A 4 5-(k4) . (6.52) 

= ^E^(P)^E5^(^,-,-^-p,-,-a)^-p,-,-^5-(k4) • (6-53) 
p k +Vf/ 

It holds that 

G*(k 1 ,k 2 ,k 3 ,k 4 )= / d P5 _(k 4 )— ^ (6.54) 

Note that the expansion of is very similar to the expansion of G\ , except for the presence of 
a special vertex associated to J. The proof of the bound (6.46) is essentially identical to the one 
for the spinless case of [BM3], to which we refer for the technical details. 

7. Correlation functions 



Once that the multiscale analysis of the partition function is completed, it is possible to apply 
the same ideas and methods to the Grassmann integrals giving the Schwinger function or the 
correlations; as the analysis is essentially identical to the one in [BM1], we will give only the main 
ideas referring to §5 of [BM1] for details. 

The density-density correlation can be written in terms of a Grassmann integral in the following 

way 

< p(x)p(y) > T -< p(x)p(y) > - < p(x) X P (y) >= ^( x )^( y ) (7-1) 

where 

5(0) = log J P(^) e - V -E„ / ^W<A",. (7.2) 

We shall evaluate S in a way which is very close to that used for the integration of the partition 
function in §2. We introduce the scale decomposition described above and we perform iteratively 
the integration of the single scale fields, starting from the field of scale 1. 
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After integrating the fields ip^\ ...tp^ h+1 ^ we find 

e s(4>) = e -Lf3E h +s<- h+1 H.<t>) J p Zh C(i ( d7 /,<ft) e -v <h) (v%:^ <£h) )+e( h »(v%:V'^ h) » ) (7.3) 

where Pz h ,o h ,c h {dip ( < h) ) and V ft are given by (3.15) and (3.16), respectively, while S^+V ((f)), 
which denotes the sum over all the terms dependent on <j) but independent of the ip field, and 
B^ h \i>^ h \<t>), which denotes the sum over all the terms containing at least one </> field and two ip 

fields, can be represented in the form, if J dx = J 2 ll dxo X^eA 

~ T 

oo „ m 

^ (/l+1) (^)=E /rfx 1 ... ( ix ro ^+ 1 )(x 1 ,...,x m )[n^(x i )] (7.4) 

to=1^ i=l 

B^\^ h \ </>) - E E E / rfx i • • • dx «> rf yi • • • ^ • 

771—1 n— 1 (7,U> ^ 

m 2n 

• < ) 2 „,Jx 1 ,.,x m ;y 1 ,... 1 y 2 „)[lI^)] [U^tt 



(7.5) 



Since the field <f> is equivalent, from the point of view of dimensional considerations, to two ip 
fields, the only terms in the r.h.s. of (7.5) which are not irrelevant are those with m = 1 and 
n = 1, which are marginal. Hence we extend the definition of the localization operator C, so 
that its action on B^ h \ip^- h \ <j>) in described in the following way, by its action on the kernels 

1) if m = 1, n = 1 then 

^^(p; k l5 k 2 ) = < 2 ^(0; 0, 0) (7.6) 

2) £ = in all the other cases 
It follows that 

CB™ (^ h) , 0) - ^-F^ h) + ^-F^ h) , (7.7) 
where and are real numbers, such that z[^ = = 1 and 

*P> = W rfx0(x)e 2 ^P-^)+V^- , (7.8) 
= E / rfx0(x)^^- . (7.9) 

<T = ±1 J 

By using the notation of §2, we can write the integral in the r.h.s. of (7.3) as 

e -Lf3t h 



e -L0th 



J Pz h _ 1 ,c h _ 1 (d^ h - 1) )- (7.10) 
• J P z / _ 1 (^(M) e -^ (h HV^V' (£h) )+e (h >(V^7V' <£h) ^) i 

where (^/Zj^tp^) is defined as in §3 and 

BW(V^fcirV ( - fc) ,0) =iBW(V^V ( - fc) ,0) • (7-11) 
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£( /l - 1 )( v /Z^/; ( -' l ~ 1) , <t>) and S^ty) arc then defined through 

e -v< h - 1 '(V2^V' < -^ 1) )+B (h " 1 HV^^^ ( - h_1) ,0)-i/3£ ; h +s( h »(0) = 

= J P Zt . 1 j-i(d# ) )e-^ ) (^ (£h, ^ W ^ KM '*' . (? ' 12) 

Of course also the new renormalization constants related to the density-density correlation func- 
tion obey to a Beta function equation of the form 

z (i) 

-^-l + z£\ i = l,2, (7.13) 



(0 



(1) (2) Z 1 

where z h and are some quantities of order g^. It turns out that Xvaih^-oo — ^tt = 1 + 0(U) 
while lim^^-oo Z\ = 1 + 0(U), with 771 = — fe?7 + 0(U 2 ) and & > is a suitable constant. The 
bounds for the expansion of the Schwingcr function or the correlation functions are done exactly 
as in §5 of [BM1]; to the first term in (1.9) or to the first two terms in (1.12) contribute only trees 
with only endpoints with scale < 0; the other trees have at least an endpoint at scale 1 so that by 
the short memory property they have a faster decay. 

8. The Hubbard model in a magnetic field 



We only sketch the analysis when there is a magnetic field as it is indeed very similar to analysis 
of the vanishing magnetic field case. 

The presence of a magnetic field destroys the SU(2) spin symmetry. The counterterms are intro- 
duced by the following definition 

t = t — ^ 8 a cos p F = [i + sign (a)h — v a (8.1) 

(7 

This means that V in the partition function (2.5) is replaced by 

,/3/2 ,/3/2 ,/3/2 

V = U dx o y2^+ + iP^ + iP+_iP^_+ / dx o y2v <T ip+ ip~ + / dx V 5 a t Xi yip+ ip~ 

_(8.2) 

The ultraviolet and infrared integration are done as in §2, §3, with the difference that for h < h only 
quartic momomials verifying || 5Z i=1 EiLOiPp 1 1| = (instead of (3.18)) are present in the effective 
potential. The definition of £ on the quartic terms is similar to (3.19) with the difference that the 
delta function in the (3.19) is replaced by S(Y^t = i ^i^iP'p )- This means that the quartic marginal 
terms verify y~]*_ 1 EiUiip'p = mod. 2ir and this condition forbids the configuration of u given by 
the second of (3.22), if h is small enough, as Pf(&) — Pf{— c) + mi 7^ 0; in other words there is no 
the analogue of the g^-terms in the effective potential. Moreover we are assuming in Theorem 2 
that I cos^ 1 (fi + h) + cos^ 1 (^, — h) — ir\ > C for some constant C; this implies that the configuration 
(u>, —u>, <j, —lu) is not allowed. 

The wave function renormalization depends by a, that is Zh.u) and the relevant part of the 
effective potential is given by 

£VW(V^/> ( -' l) ) = z Z^ F -A^Z^'i^- h) +l h Vh,*F n>a {^7 1 ^) (8.3) 

There are then in the h 7^ case 4 quadratic and 6 quartic running coupling constants; Theorem 3 is 
still valid if they are small enough. We can choose v ai S a so that Vh,+ , Vh,->8h,+ : &h- ar e 0(gh^ eh ); 
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this is shown by a fixed point argument essentially identical to Lemma 2. The four quartic running 
coupling g\ obeys to equations of the form, if i = (2o+), (2o— ), (2p+), (2p— ), (4+), (4—) 

where is given by the sum of trees with no endpoints at scale i, only g\ propagators and 
no endpoints to which are associated Vk,dk', if Uh t+ ,Uh-,8h,+,6h - are O(ghj 0h ) then, as in §5, 
K < Cg h g 9h . 

The flow of the quartic running constants is even simpler as the one in the h = case as \(3 l h \ < 
Cg~h,g 9h - This can be proved as in §6 introducing the following reference model, with 

Vl= J2 dkl - / ^E^m^^c^^^^ 

+54, < T^k 1 , w , CT C 2 , w , CT ^k 3 , W ,-aC 4 , w ,-a] ( 8 - 5 ) 

In this reference model the interaction has five independent parameters, instead of three as 
in the previous case. We can analyze by RG the reference model and we get the couplings 
9h 2 +> 9h 2 -> 9h 2 +> 9h i -> 9h 4 +> 9h 4 -• The proof that their values remains close to the initial 
value is essentially identical to the analysis in §6; (8.6) is replaced by 



< V'k 1 , + , CT ^k 2 , + ,<T , /'k 3 ,-,->k 4 ,-,- ( T >T = 

.9-,-<r(k 4 ){G : i ! _ ff (k3)[ 5 ^ (T < ^ u+itr ^k 2 ,+,aPki-k2,+,T >T 
+92,-a < V'k 1 ,+,aV'k 2 ,+, CT Pki-k2,+,-a >T 
+94,-a < V'k 1 ,+,aV'k 2 ,+, CT Pk 1 -k 2 ,-,- ( r >t] 

+ J dp[g^ a < <, + > k2 ,+,^k 3 ,->k 4 -p,-^P,+,- >t + (8.7) 

/ d P92,-a < V'k 1 ,+,aV'k 2 ,+,aV'k3,-, I TV'k 4 -p,-, CT Pp,+,-^ >T + 

J dp 34 ,-<T < i'i 1 ,+,A^,^i 3 ,+^k i -p,+^pp,+,-'r >t}} 

By using the Ward identities of §7 one gets that gX,2,+>ffi,2,-iK,2,+,9h,2-,9hA,+>9h,4- remain 
close to the initial value, and this implies that /3 l h is asymptotically vanishing. This means that in 
presence of a magnetic field one has Luttinger liquid behaviour also with a attractive interaction; 
of course this will be true only if h is non vanishing, and it is 0(h a ) for some constant a > (in 
fact h is finite if \p F — p F a \ ^ 0) . 

9. Appendix 

9.1 Ultraviolet decomposition 

It is convenient to introduce an ultraviolet cut-off N by writing 



N 



3 [ 1 ' Ar ](x,y)=^ 3 (")(x,y) (9.1) 



n=0 

where 



ff< n W) -T^E fu.v.(k)h n (k )— J—r^r (9-2) 

P kex> — ik — t cosk + t cos p F 
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with h (k ) = x(fco) and h n (k ) = \{l n+1 k ) - x{l n ko); it holds that limjv-oo g^' N ^ (x, y) 
g( u - v -)(x i y) and, for any integer K 



We define 

V^(4>)= lim log-L [ P(d^)e v ^ [1 ' N]+ V (9.4) 

We can integrate iteratively scale by scale, and after the integration of the scales N, N — 1, .., k + 1 
we get 

V^(^) = lim log-L [ P(d^)e v ^ [1 ' k] +^ (9.5) 

TV^oo TVfe J 

It is well known that can be written as sum over trees r similar to the ones in §3 (see for 
instance [GLM] for the analysis of the ultraviolet problem in the Hubbard model in any dimension) 
each of them bounded by, if m v is the number of endpoints of type U following the vertex v on t 

C"[max([/, \u\, Itfl)]^-"^- 1 ) JJ 7 -(»«- n «')("»«-i) (9. 6 ) 

V 

One can have m v — 1 only if v is a trivial vertex following the first non trivial vertices on r; 
then the terms with m v = 1 correspond to self-contractions or tadpoles; note however that no 
divergence are associate to self-contractions as g^ n ' N \x, x) is bounded uniformly in N. Consider 
then a generic tree with all the sets P v assigned; the simple expectations over the trivial vertices 
in the tree with m v — 1 before the first non trivial vertex v can be explicitly computed, giving 
^+(<n v )^-(<n v )g( n . t M)^ x y^ rest of the tree is bounded by an expression like (9.6) with 
m v > 1, so that by summing over all the scales and the trees the bound (2.17) is found. 

9.2 Spin symmetry 

Finally the symmetry property (2.18) follows from the SU{2) invariance of the Hubbard model. A 
direct way to check this property consists in expanding the truncated expectations corresponding 
to the integration of qp u - v - in terms of Feynmann graphs. The interaction can be also written in 
the following way, making more explicit the spin symmetry of the Hubbard model 

x a a' x,<7 a x,a 

(9.7) 

As usual, the Feymann graphs are obtained representing as vertices the three addends in (9.7) 
with four or two oriented half-lines, and contracting in all possible ways the half lines with con- 
sistent orientation; it is also convenient to represent the quartic term as a couple of two half-lines 
connected by a wigghly line, representing the interaction. The value of each Feynmann graph is 
obtained associating to each line a propagator g u - v (x;y) and integrating over all the coordinates; 
the contributions from graphs with four uncontracted half lines has in general the form 

j dx!...y dx 4 <^> X2 ,>+^< (T ,VF CT %,(x 1 ,..,x 4 ) (9.8) 

In order to prove that the kernel is spin- independent, that is 

W„, a - W a ,- a (9.9) 

we note that in the Feynmann graph we can identify a line of propagators g u,v (x, y) (possibly a 
point) connecting ip£ lt(T with ip^ 2 ^, and another line connecting tp^ 3 a , with a ,; on such two lines 
there are points to which are attached wiggly lines to which are attached the fields ^ „ ip~t„ip~„; 
the crucial point is that such expression does not depend from the fact that it is connected by the 



(9.3) 
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wigghly line to a a or a' line. Hence the contributions to a and W® _ a can possibly differ only 
because in one case there is a line of propagators a and in the other case — a; but the propagators 
are spin-independent hence the values of such two contributions are identical (and independent 
from a). The same argument can be repeated to prove that W% a = Wj} _ al by performing a single 
scale integration with propagator g- h (x,y). 
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